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ABSTRACT 

We present a deep K'-hand {2.12fim) imaging of 1' x 1' Subaru Super Deep Field (SSDF) taken with 
the Subaru adaptive optics (AO) system. Total integration time of 26.8 hours results in the limiting 
magnitude of K' ^ 24.7 (Scr, 0'.'2 aperture) for point sources and K' ~ 23.5 (5ct, 0'.'6 aperture) for 
galaxies, which is the deepest limit ever achieved in the K' band. The average stellar FWHM of the 
co-added image is O'.'IS. Based on the photometric measurements of detected galaxies, we obtained 
the differential galaxy number counts, for the first time, down to K' ~ 25, which is more than 0.5 
mag deeper than the previous data. We found that the number count slope dlogN/dm is about 0.15 
at 22 < i^' < 25, which is flatter than the previous data. Therefore, detected galaxies in the SSDF 
have only negligible contribution to the near-infrared extragalactic background light (EBL), and the 
discrepancy claimed so far between the diffuse EBL measurements and the estimated EBL from galaxy 
count integration has become more serious . The size distribution of detected galaxies was obtained 
down to the area size of less than 0.1 arcsec^, which is less than a half of the previous data in the 
K' band. We compared the observed size-magnitude relation with a simple pure luminosity evolution 
model allowing for intrinsic size evolution, and found that a model with no size evolution gives the 
best fit to the data. It implies that the surface brightness of galaxies at high redshift is not much 
different from that expected from the size-luminosity relation of present-day galaxies. 

Subject headings: cosmology: observations — galaxies: formation — galaxies: evolution — infrared: 
galaxies — techniques: high angular resolution 



1. INTRODUCTION 

The process of galaxy formation and evolution is 
one of the most important unsolved problems in astro- 
physics. Deep imaging of blank field is a vital method 
for studying the prope rties of g alaxies in the distant uni- 
verse (Yoshii & Takah arall98a() . The Hubble Deep Field 
(HDF) taken by the Hubble Space Telescope (HST) had 
revealed deep images o f the universe at optical wave - 
lengths ijWiinams et al.t l996. 2000; Gard ner et al.l2000j) . 
and provided us with valuable information of the distant 
universe. However, for high-redshift galaxies at z > 1, 
the rest- frame optical light, that exhibits the fundamen- 
tal structure of s tellar component i n galaxies, shifts into 
the infrared (e.g. iSheth et alJl200.'l) . Therefore, the near 
infrared (NIR) deep imaging becomes important to study 
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the formation and evolution of galaxies at high-redshifts. 
Moreover, the uncertainties due to the evolution of galax- 
ies and extinction by interstellar dust are less significant 
in the longer wavelength. Thus, deep imaging in the K 
(2.2/im) band, that is the longest wavelength at which 
the high-sensitive observations can be carried out from 
the ground, is essential to study the fundamental prop- 
erties of high-redshift galaxies. 

A number of ii'-band imaging surveys were car- 
ried out using ground-based large telescopes with 
different spatial coverages a nd limiting magnitudes 
(e.g. ICa,rdner et alJ iTflM ll^ »C 1a,zebrook ^Tal 
19941 iMcleod et all 1X9951 iDiorgovs ki et all Il99r 



Huang et al."1997^ 'Moustakas et al. 1997 : iSzokolv et alJ 
^98; Minez a ki et a l. 1998a; Bcrsha dv et alJ 119981: 
Sara cco"eraLl Il999t iV aisanen et all 120001: iMartinil 



Baker et al 



200H iMaihara et al. | 
Cristobal-Hornillos et al.l 
Among these surveys, the deep K' (2.12/im) imaging of 
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the Subaru Deep Fiel d (SDF.lMaihara et al ll200H) using 
the Subaru/CISCO (|Mntohara et a] .1120021) achieved a 
limiting magnitude oi K' ~ 23.5 ( 5a) for point sources 
with integration time of 9.7 hours. iTotani et alJ l|2001a[l 
studied the galaxy number counts in the SDF and 
suggested that a number evolution may occur or a new 
population may emerge in the faint end, while the num- 
ber count at brighter magnitude range is consistent with 
the pure luminosity evolution (PLE) without number 
evolution. 'Tot ani et al.l l|2001b(l derived the contributed 
flux of galaxies to the extragalactic background light 
(EBL) in the K band using the SDF number count 
data. They concluded, from the observed flat slope of 
differential galaxy number count {d\ogN/dm ~ 0.23) 
in the faint end and the theoretical estimate of the 
number of missed galaxies in the SDF survey, that 
more than 90% of the galaxies contributed to the EBL 
in the K band has already been resolved into discrete 
sources of galaxies. However, the EBL flux derived 
in this way accounts only for less than a half of the 
EBL flux measured in the same K band in the form 
of diffuse emission llGoriian et all [2 000: Wright 200^ 
IMatsumoto et al.ll2001t ICambresv" al...2001.) . indicat- 
ing a problem of missing K light in the universe, which 
cannot be explained by normal galaxies. Therefore, even 
deeper imaging is required to investigate the galaxy 
population in the faint end and the origin of missing 
K light in the EBL. Recently, a deep imaging of the 
Hubble Deep Field South (HDF-S) was carried out in 
the K, band (2.16 ^m) using the VLT/ISAAC (FIRES, 
^abbc ct al. 2003) and reached a limiting magnitude 
of 23.8 mag (5cr) with integration time of 35.6 hours. 
However, because unrealistically long integration time 
is required to increase the sensitivity further from this 
level, the sensitivity of deep imaging observations under 
usual seeing condition has almost reached the attainable 
limit. 

To push the limit of deep NIR imaging, we initiated a 
new deep imaging program in the K' band using the Sub- 
aru adaptive optics (AO) system. AO compensates the 
disturbed wavefront by earth's atmosphere and provides 
nearly diffraction-limited spatial resolution. Because the 
flux in the diffraction-limited core largely increases, it is 
expected to improve the sensitivity of detecting faint ob- 
jects with AO. Although the sensitivity gain with AO is 
known to be small for extended objects such as galaxies, 
we can improve the sensitivity of detecting high-redshift 
galaxies because they are expected to be compact ac- 
cording to the prediction of hierarchical cold dark matter 
(CDM) model. Studies of high resolution deep imaging 
by the HST/NICMOS have also shown that faint galax- 
ies at high redshifts are quite compact in the H (1.6/im) 
band (,Yan et al...l998.) . 

First objective of this program is to investigate the 
galaxy population in the unprecedented faint end (K' > 
24) with high detection completeness. This may reveal 
a significant faint galaxy population that contributes to 
the EBL in the K' band. Second objective is to study the 
morphology of high-redshift galaxies at rest-frame optical 
wavelengths with high-resolution AO image. Formation 
of the Hubble sequence is thought to have occurred at 1 < 
z < 2 ijKaiisawa fc Yam ada 2001). High-resolution deep 
K' imaging of galaxies in this rcdshift range may clearly 
reveal the process of morphology formation of galaxies. 



In this paper, we report the first results of our deep 
K' imaging program with the Subaru AO system. The 
layout of the paper is as follows. The field selection and 
observational strategy are summarized in §2. The pro- 
cedure of data reduction is described in §3. The proce- 
dure of source detection and photometry is described in 
§4. The analysis of reduced images and detected source 
counts is given in §5. The quality of our AO deep K' 
imaging is examined in §6. The results of the galaxy 
number counts, the contributed fiux of galaxies to the 
EBL, and the size distribution of detected galaxies are 
discussed in §7. Finally, the summary is given in §8. 
Throughout this paper, we adopt the magnitude system 
where Vega is 0.0 mag, and the cosmological parameters 
of f7M=0.3, f}A=0.7, and iJo=70 Km/s/Mpc. 

2. OBSERVATIONS 

Our observed field is a part of the well-studied deep 
field called "Subaru Deep Field (SDF)", which is a blank 
sky r egion near the north galactic pole (*M aihara et alJ 
1200 1|) . This field h as been observe d extensively by the 
Subaru Telescope l|Ive et all 12004(1 as an observatory 
proje ct with a wide wav elength coverage fr om op tical to 
NIR ((Maihara et alJl20?) l: Kahikawa ct al] l200l . JK'- 
band images for a 2' x 2' area o f the SDF were obtained 
with the Subaru/CISCO (Motohara et al. 2002) and 
BV Ri' z' -hand images for a 30' x 37' area, which includes 
the CISC O field, were obt ained by the Subaru/Suprime- 
cam ( Miv azaki et al.l2002|) . The SDF was originally cho- 
sen such that a reasonably bright star is located in the 
field as a wavefront guide star for AO observation. We 
placed this bright (i? ~ 12) star near the center of our 
1' X 1' AO field. The center coordinates of our AO 
field was set to be a=13''24"23!6 and (5=+27°30'30'.'4 
(J2000). We call this field "Subaru Super Deep Field 
(SSDF)". The location of the SSDF relative to the SDF 
is shown in Figure^ Because very deep BVRi'z' images 
of the SDF have already been obtained, it is possible to 
determine the photometric redshift of detected galaxies 
in the SSDF. Seeing-limited infrared J and K' images of 
the SDF have also been taken by Maihara ct al. (200ll) 
and a part of this field is overlapped with our AO field 
(see Figure^. Thus, the quality of our AO deep im- 
ages can be directly compared to that of conventional 
seeing-limited deep images (sec §6). 

Our observ ations were carried out using the Subaru 
AO system (Takami et al.ll2004() and the IRCS (Infra- 
Red Camera and Spe ctrograph, iTokunaga et al.l Il998t 
iKobavashi et al.lf2000() both mounted on the Cassegrain 
focus of the Subaru Telescope. The IRCS imager is 
equipped with a Raytheon Aladdin HI 1024x1024 InSb 
array, offering two pixel scales of 23 and 58 mas for 
imaging in a wavelength range of 0.9—5.5 /xm. We used 
the pixel scale of 58 mas to cover a wider field of view 
with the K' filter (1.96-2.30 /im, Fx = 4.66 x 10"!° 
Wm-2/ini-\ 'Tokunagact al. 2002). The Subaru AO 
system uses a curvature sensor with 36 control elements, 
which provides a stellar image of Strehl ratio (the ratio 
of the observed star peak to the peak value of a perfect 
telescope diffraction pattern) of '^0.28 and full-width at 
half maximum (FWHM) of ~ 0'.'07 in the K band with 
an i? ~ 12 guide star under best observing condition. 
The improved image quality is expected to increase the 
detection sensitivity to a point source by 1 mag. How- 
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Fig. 1.— The location of the Subaru Super Deep Field (SSDF) 
relative to the Subaru Deep Field (SDF). Intensive observations 
have been carried out in the SDF using the Subaru Suprime-Cam 
(BVRi'z') and the CISCO (JK'). The stellate symbol in the SSDF 
shows the location of the AO guide star (R ~ 12). 



TABLE 1 
Summary of K'-band observations. 



Date (UT) 


Exp. time [sec] 


Field 


2003/03/17 


18630 


SSDF 


2003/03/18 


18000 


SSDF 


2003/03/19 


10800 


SSDF 


2003/03/20 


18480 


SSDF 


2003/04/22 


18630 


SSDF 


2003/04/23 


18630 


SSDF 




135 


M13 (PSF reference) 


2003/04/24 


8100 


SSDF 




135 


M13 (PSF reference) 



ever, because AO correction performance degrades with 
increasing distance from a guide star (see §6.3 for de- 
tails), the Strehl ratio at the area away from a guide 
star is lower than the best value. Moreover, since the 
Strehl ratio should vary with time depending on the ob- 
servational condition, the resultant Strehl ratio for long 
exposure time is averaged and lower than the best value. 
Similarly, the FWHM for long-time exposure image could 
be broader than the best value at the area away from a 
guide star. For the present data, the achieved Strehl ratio 
and the FWHM, which were measured from a point-like 
source located at 24 arcsec away from the guide star, 
are 0.1^ and O'.'IS on average, respectively. Because the 
galaxies in our field of view are distributed within about 
40 arcsec away from the guide star, our measured Strehl 
ratio and FWHM are typical for our observation. 

We carried out our observations for a total of nine 
nights with almost photometric condition. The observ- 
ing log is given in Tabled During the observations, we 
repeated the set of nine-position dithering with a 3x3 
grid of 2'.'0 separation. Individual exposure time of 90 or 
120 sec was adopted such that the sky background does 
not saturate the detector well. Since a few bright stars 
reached saturation even with this short exposure time, 

^ Since the pixel scale of 58 mas is not small enough to sample 
the small size of point spread function, the estimated Strehl ratio 
may not be accurate. 




Fig. 2. — The variation of the sky background brightness in the 
K' band (a) and the spatial resolution with AO (b) during the 
observations. Different symbols correspond to different observing 
nights listed in Table [Tl The spatial resolution is represented by 
the FWHM of the point-like source ("S" in Figure |2J. 

residual images following the dithering pattern appeared 
around the stars due to the memory effect of the InSb 
array (see grid patterns around bright stars in Figure |3J). 

Throughout each observing night, the variation of the 
sky background brightness was small except at large air- 
mass. The sky variation corresponds to about ±0.1 mag 
in the limiting magnitude variation of each frame (Figure 
suggesting that the observing condition was almost 
stable throughout all observing nights. The correction 
performance of AO changes every moment depending 
largely on the observational condition. We monitored the 
variation of spatial resolution by measuring the FWHM 
of a point-like source which has the sharpest and nearly 
circular profile in our field of view (marked "S" in Figure 
13) except for the saturated bright stars. Since this point- 
like source "S" is too faint to be detected in each single 
frame, we combined each set of nine-point dithering to 
measure the FWHM. The variation of FWHM after AO 
correction was found to be small for all frames (Figure 
Eb)- Thus, the correction performance of AO was sta- 
ble and the point spread function (PSF) remained to be 
constant for each night. 

3. DATA REDUCTION 

We reduced the data with IRAF^ software packages. 
Before the data reduction, we checked all data frames vi- 
sually and removed some low-signal frames due to cirrus 
as well as some frames with unusual dark patterns due to 
unstable detector temperature at the beginning of contin- 
uous exposures. We rejected these frames (~ 15% of all 
frames) before the data reduction and the resultant total 
integration time of all used frames is about 26.8 hours. 
Flat fielding was performed using the median combined 
sky fiat frame. The sky flat frames were created for each 
night using the dark subtracted raw frames. In generat- 

IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Fig. 3. — IRCS+AO /sT'-band image of the Subaru Super Deep Field (SSDF) with the total integration time of 26.8 hrs. The field size is 
1' X 1' with the pixel scale of 58 mas. The stellar FWHM was measured to be 0'.'18 for the point-like source (indicated by "S" mark). The 
bright star at the center of the image was used as the AO guide star. Ghosts of bright stars due to the IRCS internal optics (beam splitter 
and compensator) are indicated by cross marks. The grid patterns around the guide star and other bright objects are the residual image of 
previous frames due to saturation of the detector. Because of these mock objects, some areas near the bright stars and their ghosts (boxed 
areas) were excluded from the detection area. Gl is a disk galaxy with an effective radius of 0'.'25, t hat is the average size of galaxies in 
the SSDF image (see §6.2). G2 is a galaxy whose radial profile and model fitting are shown in Figure [T^ 



ing the sky flat frames, bright and extended objects in 
each raw frame were masked in order to reduce their con- 
tribution to the final flat frame. We made the position 
map of the bad pixels on the detector from the dark and 
flat frames, which were used for picking the pixels with 
extraordinarily high counts and the pixels with extremely 
low sensitivity, respectively. Bad pixels on our observed 
frames were removed by interpolation using this position 
map. The nine frames of each dithering set were median- 
combined to create the sky frame for each set. Before this 
process, all of the discernible objects are removed with 
interpolation by adjacent sky counts to reduce the contri- 
bution of the individual objects to the resultant sky level. 
The sky subtraction was performed for each set of nine- 
point dithering (every 810 or 1080 seconds) to minimize 
the effect of time variation of the sky background. Then, 
we shifted and combined all of the sky subtracted frames 
with exposure time weighting. The image offsets were 



determined at a sub-pixel level from the bright objects 
in each frame. The final K' image is shown in Figure 
13 The average FWHM of final image was about O'.'IS 
for the relatively bright [K' ~ 21.2) point-like source 
"S" . Throughout this paper, we used the profile of this 
point-like source as the PSF of our AO image. 

4. SOURCE DETECTION AND PHOTOMETRY 

Source detection and photometry were performed by 
the SExtractor (Bertin & Arnouts 1996). Before the 
source detection, the final image was smoothed by a gaus- 
sian filter (a = 1 pix) to give an optimal source detection 
capability with less spurious detection. We defined the 
detection threshold as the 1.5ct level of the surface bright- 
ness fluctuation of the sky (23.64 mag/arcsec^). If more 
than 19 contiguous pixels have larger counts than the 
threshold, we regarded it as a positively detected source. 
To reduce the risk of spurious detection, some areas near 
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Fig. 4. — Measurement of background noise within the aperture, (a) Gaussian function (dotted hne) is fitted to the histogram of measured 
flux in the aperture area A = wr^ (r = 10 pixels) which was randomly placed avoiding the object locations. We regarded the width cr 
of gaussian function as the background noise within the aperture. This method allows us to measure the true background noise in which 
a pixel-to-pixel correlation is taken into account, (b) Relation of aperture radius r versus background noise within the aperture. Filled 
circles show the measured background noise and the solid line shows the fitted line to the data with a two dimensional function, ar^ + br, 
where a and b are the free parameters. The measured background noise with large aperture is significantly larger than that expected from 
a linear scaling of the pixel-to-pixel noise (dashed line), probably due to the correlated fluctuations of the background on large scale. 




Aperture Radius [arcsec] 



Fig. 5.— Signal-to-noise ratio (SNR) of the K' ~ 21.2 point- 
like source "S" (solid line with fllled circles, left y-axis) and the 
K' ~ 20.8 disk galaxy "Gl" (dashed line with open circles, right 
y-axis) as a function of aperture radius. This diagram is used to 
determine the 5<j limiting magnitude. Dotted line shows the maxi- 
mum SNR and the aperture radius. In case of the point-like source 
"S", maximum signal-to-noise ratio of 127 is achieved at aperture 
radius of O'.'l, so that the 5(t limiting magnitude for point sources 
is 21.2 - 2.51og(5/127) ~ 24.7 with an aperture diameter of 0'.'2. 
Similarly, the 5cr limiting magnitude for galaxies was estimated at 
K' ~ 23.5 with an aperture diameter of 0'.'6. 

the bright stars and the known ghosts of bright stars are 
excluded from the detection area (see Figure |3J). As a re- 
sult, the total number of detected sources is 236 within 
the net detection area of 0.9 arcmin^, although some spu- 
rious detections due to statistical noise may be included 
in this number. 
We performed photometry of detected sources in terms 



of aperture magni tude, isop hotal magnitude, and total 
magnitude (see iBertin fc Arn outs (1996) for further de- 
scription of each magnitude). Throughout this paper, 
we mostly use the total magnitude except for the dis- 
cussion of size distribution in §7.2, where we use the 
isophotal magnitude. The photometric calibration was 
carried out using the in frared faint standard star GSPC 
P330-E persson et aU lTQOS). We took the image of this 
standard star at airmass ~ 1 in a photometric night and 
used it as the photometric reference for the object frames 
which were taken in the same night at the same air- 
mass. Then, a bright galaxy {K' ~ 16.9) in the field 
was used as the photometric reference for the final im- 
age. The resultant zero point in the final combined image 
is K' - 23.645 ± 0.011 mag/ADU. 

Understanding the noise properties is crucial because 
the limiting magnitude and photometric uncertainty rely 
on them. The photometric uncertainty in typical near- 
infrared deep images is well described by the poisson 
noise of the background signal in each pixel. However, 
image processing, such as shift and combination, has 
introduced correlations between neighboring pixels. If 
the photometric uncertainty was estimated from linear 
scaling of the pixel-to-pixel variation, we would under- 
estimate it. To estimate the true uncertainty of pho- 
tometry, we measured the sky fluxes in several tens of 
circular apertures at random position in the final im- 
age and regarded the standard deviation a of the mea- 
sured fluxes as the photometric uncertainty. We changed 
the radius of circular apertures ranging from 1 to 40 
pixels (0'.'058— 2'.'3) to cover all the photometric aper- 
tures of detected sources (Figure 0}. We fitted a with a 
two-dimensional function of the aperture radius, cr(r) = 
ar^ + br, where a and b are the fitting parameters. We 
used this function for calculating the photometric uncer- 
tainty of each detected sources with various apertures. 
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The measured uncertainty with this method largely ex- 
ceeds the uncertainty expected from linear scaling of 
the pixel-to-pixel noise (dashed line in Figure ^ . The 
large discrepancy a t large aperture was also reported in 
iLabbe et alJ (|2003) . This could be caused by correlated 
fluctuations of the background on large spatial scale. 

We estimated the limiting magnitude for point sources 
using the point-like source "S" in the final image. The 
limiting magnitude was determined by using a diagram 
that shows signal-to-noise ratio (SNR) of the point-like- 
source "S" as a function of an aperture radius (solid 
line in Figure EJ, that was represented as a ratio of 
flux to photometric uncertainty within an aperture. The 
achieved 5a limiting magnitude for point sources is K' ~ 
24.7 (5(t) with an aperture diameter of 0'.'2, where the 
aperture size was determined to have the maximum SNR. 
This is the faintest i^-band limiting magnitude for point 
sources achieved to date. Because galaxies would have 
more extended profile than point sources, the limiting 
magnitude for galaxies should be brighter than that for 
point sources. Thus, we also estimated a limiting mag- 
nitude for galaxies with the similar method as for point 
sources using a galaxy "Gl" (see Figure|3l, that is a disk 
galaxy which has the typical effective radius {r^ — 0'.'25) 
in the SSDF (see §6.2 for details). The achieved 5cr lim- 
iting magnitude for galaxies is K' ^ 23.5 (5ct) for 0'.'6 
aperture. The number of detected sources brighter than 
the limiting magnitude for point sources is 145, while the 
total number of detected sources is 236 with the faintest 
magnitude of K' 26.0. 

5. ANALYSIS 

5.1. Detection completeness 

At the fainter magnitude, the detection completeness 
decreases because of the statistical noise. We estimated 
the detection completeness for point sources by conduct- 
ing the same source detection as described in §4 but with 
a large number of artificial point sources placed at ran- 
dom positions in the final image. The artificial point 
sources were created from the point-like source in the fi- 
nal image (source "S" in Figure by applying a flux 
scaling in the range of 20 < K' < 27. Thick solid line 
with flUed circles in Figure El shows the detection com- 
pleteness curve for point sources in our flnal image as 
a function of K' magnitude. We estimated that 50% 
completeness for point sources is achieved at K' ~ 25.0, 
which means a half of all galaxies can be detected at 
this magnitude, and it is about 0. 6 magnitude deeper 
than previous deep imaging such as ILabbe et al.l l|2003fl 
and ^aihara et al. (2001). Thus, our data should of- 
fer the most reliable source detection down to 25 mag. 
We also estimated the detection completeness for galax- 
ies with artiflcial galaxies that have the typical size in 
the SSDF. The artiflcial galaxies were created to have 
the same surface brightness profile as the galaxy "Gl", 
which was used in the limiting magnitude estimation in 
§4. We placed them at random positions with varied flux 
in the range of 20 < K' < 27. Figure [3 shows the es- 
timated detection completeness curve for galaxies. The 
completeness curve for galaxies in our final image are 
shown as thick solid line with filled circles. For a compar- 
ison, the completeness cur ve for extended sourc e in the 
previous deep imaging by iMaihara et"all l|200H) is also 




K' total magnitude 



Fig. 6. — Detection completeness curves derived from the simu- 
lations using the artificial point sources. Thick solid line with filled 
circles shows the detection completeness for point source of our K' 
survey in the SSDF, while thin lines with other symbols sh ow that 
of deep K or K' imaging surveys from the literature I Labb^^^al 
2003; Maihara et al.li 2001t IBershadv et al.l[T998i: .Moustakas et al 
1997; Diorgovski ct al? lS9a) . Horizontal dashed line indicate the 
level of 50% completeness. 

shown as thin dashed line with stellate symbols in the 
same diagram. We estimated that 50% completeness for 
extended source is achieved at K' c:^ 24.1, that is about 
0.6 mag fainter than that of M aihara et alJ l)2001|) . This 
magnitude gain is same as for the gain for point source, 
suggesting that the AO system improves the efficiency of 
detecting high-redshift galaxies as well as that of point 
sources. In this paper, the detection completeness for 
point source is used to correct the number counts. There- 
fore, the detection completeness here should be consid- 
ered upper limits and then the corrected counts should 
be considered lower limits. 

5.2. Noise contamination to number counts 

We detected the objects down to K' < 26 in our fl- 
nal image. Based on the photometric measurements of 
these objects, we derived the number counts in the 0.5 
magnitude bin (Table The number counts, in princi- 
ple, contain not only the counts from galaxies, but also 
the counts from the Galactic stars and spurious detec- 
tions due to the effect of statistical noise. The contribu- 
tion of the Galactic stars to our number counts should 
be negligible, because the expected star counts toward 
the Galactic pole are less than 0.2 mag~^ acrmin~^ in all 
magn itude range ijMinezaki et al. 1998b; F ueal fc Moodvl 
12003(1 . On the other hand, the contribution of spurious 
detection is signiflcant in the faint end, where the signal- 
to-noise ratio for the detection becomes lower. We es- 
timated the number of spurious detections from an ar- 
tificial noise frame which was created by combining the 
sky subtracted frames without adjusting the dithering 
offset. The combination was conducted with the options 
of exposure time weighting, which is similar to the com- 
bination of the final image. Before the combination, all 
discernible objects in each image were removed by inter- 
polation and the estimated noise from the adjacent pixels 
was added on the interpolated pixels. Even after remov- 
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Fig. 7. — Detection completeness curves derived from the sim- 
ulations using the artificial galaxies that have the typical size in 
the SSDF. Thick solid line with filled circles shows the detection 
completeness for galaxies of our K' survey in the SSDF, while 
thin d ashed line with stellate symbols show that ofjvlaihara ot al. 
^2()(!i^ for a comparison. Horizontal dashed line indicate the level 
of 50% completeness. 



_ 100 p— ^ 1 , 1 , 1 , 1 ^ 

ft» " 




a 

I I I I 1 I I I I L_ 

20 21 22 23 24 25 
K' total magnitude 



Fig. 8. — Galaxy counts in the SSDF in the 0.5 magnitude 
bin without any correction except for subtraction of the noise con- 
tamination. Error bars show the uncertainties due to the poisson 
statistics. 



ing the discernible objects, even fainter objects could be 
remained in the noise in each single frame. In order to 
minimize the contribution of the light from such faint 
objects, the sign of each image was reversed before the 
combination. We applied the same detection procedure 
to the artificial noise image as used for the source de- 
tection. The resultant number of spurious detections in 
each 0.5 magnitude bin is given in Tabled We derived 
the galaxy number count in each magnitude bin by sub- 
tracting the number of spurious detection from the raw 
count. The estimated galaxy counts are listed in Table |21 
and shown in Figure |S| 



5.3. The scatter of number counts due to photometric 
uncertainty 

Even if the noise counts are removed, the scatter of 
galaxy number counts among magnitude bins still re- 
mains due to photometric uncertainty. In order to eval- 
uate the effect of this scatter, we performed a simula- 
tion by adding artificial objects to the final image and 
applying the same detection procedure as described in 
§4. In this simulation, about 10 artificial objects were 
added 450 times (total of about 4,500 objects) in or- 
der not to largely change the number density of the 
objects in the final image and the detection procedure 
was performed each time. The artificial objects of var- 
ious magnitudes were created from the point source in 
the final image by scaling its flux. The magnitudes of 
the artificial objects were assumed to distribute with 
power-law, N{m) cx 10"™, where N is the number of 
the artificial objects in each 0.5 magnitude bin, m is 
the central magnitude of each magnitude bin, and a is 
the power-law index. We performed the simulation only 
with point sources because previous analyses related to 
the number counts are based on point sources to give 
the lower limit of the counts. An initial power-law in- 
dex of ao — d log N/ dm was assumed to be the slope 
of raw galaxy counts over a reliable magnitude range, 
22 < K' < 24, where the signal-to-noise ratio is greater 
than 3 (i.e. more than 10 raw counts) and the detection 
completeness for point sources is higher than 90%. 

Based on the simulation, we first generated the transfer 
matrix Ty, each element of which gives the fraction of 
galaxies with magnitude mj but detected with m^. Then, 
we generated the probability matrix Pji, each element 
of which gives the probability that an artificial object 
detected with actually has a magnitude of mj . These 
elements are described as 

Pji = TijUj/ '^TikUk , (1) 

k 

where Uj is the number of the artificial objects with mag- 
nitude nij . The probability matrix is shown in Figure El 
With this probability matrix, we corrected the galaxy 
number counts as 

n- = ^P,,nf, (2) 

i 

where n™'' is the corrected galaxy number count at mag- 
nitude ruj , and nf ' is the raw galaxy count with detected 
magnitude mi (Table |2J|. If the slope of corrected galaxy 
counts was different from the initial slope ao used in the 
simulation, the correction procedure described above was 
repeated with different initial slope until the corrected 
slope coincides with the initial slope. This method has 
been employ ed in some studies of galaxy number counts 
llSmai] et aLIII 99^ IMineza,ki et al.111 998aHMa,ihara et alJ 
12001 ! 

We derived the differential galaxy number counts 
(7Vf°'') in the unit of number mag~^ deg~^ from the 
counts n'^"^ after correcting for the detection complete- 
ness estimated for point source in §5.1 (Tabled. Fig- 
ure^! shows the plot of galaxy number counts with and 
without the completeness correction. At brighter mag- 
nitude [K' < 22.0), our differential number counts are 
not complete with less than 3cr measurements (i.e. less 
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TABLE 2 

The corrected /C'-band galaxy number counts. 



K' 






rigal 






Completeness 


^cor 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


20.25 


2 





2 


1.41 


1.51 


0.96 


1.959E4 


1.287E4 


20.75 


5 





5 


2.24 


2.59 


0.99 


4.837E4 


2.157E4 


21.25 


13 





13 


3.61 


4.97 


0.98 


9.735E4 


4.168E4 


21.75 


8 





8 


2.83 


3.54 


0.95 


9.112E4 


3.067E4 


22.25 


18 





18 


4.24 


6.39 


0.96 


1.360E5 


5.479E4 


22.75 


11 





11 


3.32 


4.44 


0.94 


1.253E5 


3.871E4 


23.25 


22 





22 


4.69 


7.47 


0.96 


1.598E5 


6.367E4 


23.75 


17 


1 


16 


4.24 


6.06 


0.92 


2.100E5 


5.377E4 


24.25 


37 


4 


33 


6.40 


11.15 


0.88 


2.497E5 


1.038E5 


24.75 


31 


8 


23 


6.24 


9.47 


0.70 


2.577E5 


1.104E5 


25.25 


42 


29 


13 


8.43 


10.30 


0.33 


3.038E5 


2.587E5 


25.75 


21 


13 


8 


5.83 


8.32 


0.12 


4.746E5 


5.575E5 



Note. — (1) K -band total magnitude. (2) Raw counts of detected objects in the 0.5 
magnitude bin. (3) Noise counts in the 0.5 magnitude bin. (4) Galaxy counts in the 0.5 
magnitude bin {ugai — riraw — nnoise)- (5) Uncertainties in the galaxy counts coming from 
poisson statistics. (6) Uncertainties in the galaxy counts coming from poisson statistics 
and clustering of galaxies. (7) Detection completeness for point source. (8) Differential 
galaxy number counts in mag~^deg~^ corrected for the incompleteness and the scatter due 
to the photometric uncertainties. (9) Uncertainty in the corrected differential number counts 
coming from the poisson statistics and the clustering of galaxies. 




21 22 23 24 25 
Model magnitude nrij 

Fig. 9. — The probability matrix for detected galaxies in the 
SSDF based on the simulation using the artificial point sources. 
Each element of the matrix gives the probability that a galaxy 
detected with actually has a model magnitude of rrij. The 
probability is color-coded according to the scale bar on the right 
of this figure. 

than 9 raw counts) because of a small field of view (I'x 
1'). At fainter magnitude, the detection completeness for 
point source becomes lower than 50% at K' > 25. Thus, 
the reliable magnitude range of our differential number 
counts is 22 < K' <25. 

5.4. Field-to-field variations of galaxy number counts 

Clustering of galaxies could lead to the systematic un- 
certainties in the number counts. We estimated the field- 
to-field variations of the counts due to the clustering from 
an angular correlation function. Considering an angular 
area of fl, with a mean count of (N) galaxies, the vari- 
ance in the number of detected galaxies is 



(N) 



{Nl 




31 22 23 24 25 
K' total magnitude 



26 



(3) 



Fig. 10. — Differential ii''-band number counts of galaxies esti- 
mated from the galaxy counts in the SSDF. Filled circles show the 
galaxy counts in a reliable magnitude range of 22 < K' < 25, where 
the signal-to-noise ratio for raw counts is greater than 3 and the 
detection completeness for point source is higher than 50%. Cir- 
cled dots show the less reliable galaxy counts. These galaxy counts 
are corrected for the detection completeness and the scatter due to 
the photometric uncertainty. Open circles show the galaxy counts 
without completeness correction. The completeness correction be- 
comes significant at K' > 25. Error bars show the uncertainties 
due to the poisson statistics. 



where lu{0) is the angular correlation function of galax- 
ies and 6112 is the ai igle between the solid- angle ele- 
ments drii and df22 l|Croth fc Peebleal IToTTI) . In this 
formula, first term comes from poisson statistics and sec- 
ond term comes from the clustering. The angular corre- 
lation function for faint galaxies has been me asured in 
the K band in a number of li teratures (e .g., Baugh et alJ 
119961 iCarlberg et al., 1997: .Roche et"al][l999) , and it is 
described in a form of uj{9) = A9~^-^. The amplitude 
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Fig. 11. — The galaxy counts in the SSDF with the uncertainties 
due to the poisson statistics and the clustering of galaxies estimated 
from the angular correlation function for the 0.9 arcmin'^ area of 
the SSDF (filled circles). The galaxy counts for four discrete 0.9 
arcmin^ areas in the SDF are also plotted (open circles). The 
estimated uncertainties of the SSDF due to a poisson statistics and 
a clustering are almost consistent with the variance of the galaxy 
counts for four areas in the SDF. 

A at fixed is a decreasing function of K magnitude, 
while there is an evidence that the amplitude becomes 
relatively flat at K > 20 with a value of A ^ 1. 1 x 10~^ 
when 9 is measured in degree l|Roche et al.ll999|) . In this 
paper, we assume the amplitude-magnitude relation does 
not turn over at K > 20, which is theoretically reason- 
able (e.g., Roche & Eales 1999), to set the upper limit 
on the uncertainty of the counts coming from the clus- 
tering of galaxies. For the 0.9 arcmin^ area of the SSDF, 
the variance becomes (N) +61.3A{N)'^. The estimated 
uncertainties of the galaxy counts coming from poisson 
statistics and clustering at each magnitude bin are sum- 
marized in Table [2| The uncertainties due to poisson 
statistics and clustering are not much larger than that 
only due to poisson statistics. 

In order to check the validity of the variation estimated 
from the angular correlation function, we also estimated 
the fi eld-to-field variations of the counts using the SDF 
data ijMaihara et al.ll2001[) . which has almost four times 
larger area (1.'97 x 1.'9) than that of the SSDF. We de- 
rived the counts from four discrete 0.9 arcmin^ areas in 
the SDF at 20 < K' < 23 where the detection com- 
pleteness is higher than 90% for both of the SDF and 
the SSDF. The galaxy counts of the SSDF are plotted in 
Figure lTTl with the uncertainties estimated from equation 
(O . The galaxy counts for four discrete areas in the SDF 
are also plotted. It is found that the uncertainties due 
to the poisson statistics and the clustering of galaxies 
are almost consistent with the variations of the galaxy 
counts from four areas in the SDF. 

6. PERFORMANCE OF AO DEEP IMAGING 

6.1. Comparison with seeing-limited observation 

To estimate the performance of our AO deep imag- 
ing, we compared the quality of our AO image of the 
SSDF to that of the seeing-limited image of the SDF 



Fig. 12. — Encircled flux of the point-like source "S" in the 
SSDF (filled circles) as a function of distance from the center of 
the stellar image. For a com parison, enc ircled fiux of a point source 
in the SDF (stars. Maihara et al.ll2001l) is also shown. Since there 
is no point source in the overlapped region between the SSDF and 
the SDF, a different point source is used for the SDF. Total flux 
was normalized to unity for both point sources. Vertical error bars 
show the photometric uncertainties in each aperture. Because the 
point-like source "S" in the SSDF is much fainter than the point 
source used for the SDF, the photometric uncertainties for the 
SSDF data are much larger than that for the SDF data. Much 
higher concentration of flux was attained in the SSDF compared 
to the SDF. 

(jMaihara e t al."2001^ The SDF data were obtained by 
the Subaru/CISCO under the FWHM~0'.'45 seeing con- 
dition with the total integration time of 9.7 hours. Be- 
cause a part of the SDF field of view is included in the 
SSDF (see Figure 0, we can make a direct comparison 
between AO and seeing-limited data. 

Figure El shows the encircled flux plot of the point 
sources detected in the SSDF and the SDF. The 50% 
of the total flux is contained within a radius of O'.'IS for 
the SSDF, while 0'.'27 for the SDF. The higher flux con- 
centration due to the AO correction resulted in much 
higher sensitivity. To estimate the sensitivity gain for 
our AO observation against seeing-limited observation, 
we derived the 5(t limiting magnitude for point sources 
in the SSDF and in the SDF with the same method as 
described in §4. The 5(T limiting magnitude is about 22.9 
(0'.'2 aperture) for the SSDF and 22.4 (0'.'4 aperture) for 
the SDF with an hour integration time. Thus, the sensi- 
tivity gain due to AO correction is about 0.5 mag in the 
K' band. 

The improvement of spatial resolution with AO is 
also significant. Figure El shows a model fit to the 
surface brightness profile oi a, K' ^ 21.3 galaxy at 
z ~ 1.46 ("G2" in Figure O, using a PSF-convolved de 
Vaucouleurs profile (dc Vaucoulcurs 1948) for a bulge- 
dominated galaxy and an exponential profile l)Freemari 
[T97Q|) for a disk-dominated galaxy. While it was diffi- 
cult to distinguish these model profiles for the seeing- 
limited SDF data, this galaxy is evidently a disk dom- 
inated galaxy for our SSDF data. Effective radii of lo- 
cal galaxies, except for compact dwarf galaxies, range 
from '- ^l.O kpc to ~ 10 kpc depending on their lumi- 
nosity l|Bender et al.lll992l:llm3ev et alJll996j) . bmce our 
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Fig. 13.— Radial profile fitting to a /C' ~ 21.3 galaxy at z ~ 1.46 ("G2" in FigureEJ in (a) the SSDF with IRCS+AO and (b) the SDF 
with CISCO IJVLaihara et alj|2001l) . Filled circles show the observed data in both figures. The PSF-convolved de Vaucouleurs (dotted line) 
and exponential (dashed line) profiles were used to model a galaxy dominated by bulge and disk, respectively. Solid lines show the PSF 
profile in both figures. The high-resolution (FWHM = O'.'IS) image of the SSDF clearly shows that the best fitted model is the exponential 
(disk dominated) profile, while it is difficult to distinguish the two models with the seeing-limited image of the SDF because of the large 
FWHM (= 0'.'45). 

spatial resolution of O'.'IS corresponds to about 1.4 kpc at 
z ^ 1, we can determine the size as well as the morpho- 
logical type even at z > 1. Thus, our AO morphology 
data should enable systematic and quantitative study of 
rest-frame optical morphology of galaxies at z > 1 (Mi- 
nowa et al. 2005, in preparation). 

6.2. Effect of partially corrected PSF 

Typical AO system operated with a guide star of mod- 
erate brightness can only partially correct for turbulence- 
induced distortions. The partially corrected PSF consists 
of two components: a diffraction-limited core superim- 
posed on a seeing halo. As the Strehl ratio becomes 
lower, which means a decrease of the degree of correction, 
a greater proportion of energy in the diffraction-limited 
core is scattered into the surrounding halo and the peak 
of the diffraction-limited core becomes lower. As a result, 
the seeing halo becomes the dominant component in the 
PSF. For the present SSDF data with the Strehl ratio 
of about 0.1, since about 90% of the energy is scattered 
into the surrounding halo, the FWHM of the observed 
PSF (~ O'.'IS) is broader than the diffraction-limited core 
(~ 0'.'057 in the X-band) and the wings extended out to 
the seeing size are seen around the edge of the PSF (see 
Figure ini). 

The surface brightness profile of galaxies detected in 
the SSDF may be distorted by the extended wings of 
the observed PSF. To estimate the effect of the extended 
wings of the observed PSF on the flux and size measure- 
ment of the galaxies in the SSDF image, we conducted 
a simulation using artificial galaxies that were convolved 
with a model PSF profile. Figure E| shows the radial 
profile of the observed PSF and the fitted lines with a 
moffat profile and a gaussian profile. These profiles are 
described as 




0.2 0.4 0.6 

Distance fronn Center [arcsec] 



Fig. 14.— Radial profile of the observed PSF in the SSDF. The 
peak surface brightness of the PSF is normalized to unity. The 
moffat profile (solid line) and the gaussian profile (dashed line) are 
fitted to the observed PSF. The moffat profile represents well the 
observed PSF, while the gaussian profile only represents the central 
part of the observed PSF and failed to represents the extended 
wings. The observed PSF has a broader FWHM than the Subaru 
diffraction-limited PSF (dotted line) because of the low Strehl ratio 
(~ 0.1). 



I{r) =\l + (2^/'^ - 1) {2r/FWHMY 



(moffatlA) 



/(r)=exp -\n{2){2r/FWHMf (gaussian), (5) 

where r is the distance from the center of the PSF, 
FWHM is the full width at half maximum of the PSF 
and /3 is the constant related to the shape of the pro- 
file. As /3 decreases, the moffat profile deviates from the 
gaussian profile and comes to show more extended wings. 
For the observed PSF in the SSDF, best fit parameters 
of FWHM and (3 arc O'.'IS and 1.73, respectively We 
found that the moffat profile represents well the observed 



Subaru Super Deep Field with AO I 



11 




Fig. 15. — The isophotal area versus isophotal magnitude dia- 
gram for the observed galaxies (crosses) . We also plotted the mean 
isophotal area for the simulated disk galaxies with effective radius 
(re) of 0''.1 (open triangles), 0'.'25 (filled circles), and O'.'S (open 
squares). The horizontal error bars show the bin size around the 
mean and the vertical error bars show the Icr dispersion of the 
isophotal area of the simulated galaxies. The simulated disk galax- 
ies are convolved with the observed PSF modeled by the moffat pro- 
file. The isophotal areas of the observed galaxies are roughly dis- 
tributed between that of the simulated disk galaxies with = 0''1 
and 0"5. Thus, we used the disk galaxy with the intermediate 
size of re = 0'.'25 as the artificial galaxy that models the typical 
galaxies in the SSDF image. 

PSF, while the gaussian profile fails to represents the ex- 
tended wings of the observed PSF (see Figure^)). Thus, 
we used the fitted moffat profile to convolve the artifi- 
cial galaxies in the simulation. The artificial galaxies are 
created to have the typical size in the SSDF image. Fig- 
ure El shows an isophotal area of the observed galaxies, 
that means the area where the surface brightness of the 
galaxies is brighter than the detection threshold of 23.64 
mag/arcsec^, as a function of the isophotal magnitude. 
We also plotted the isophotal areas of the simulated disk 
galaxies with re = O'.'l, 0'.'25, and 0'.'5. Because the sizes 
of the observed galaxies are roughly distributed between 
the size distribution of the simulated disk galaxies with 
C/.'l and O'.'S, we used the disk galaxy with the intermedi- 
ate size of Te = 0'.'25 as the artificial galaxy that models 
the typical galaxies in the SSDF image. The size and 
morphological type of the artificial galaxies are same as 
the galaxy "Gl" , which was used to estimate the limiting 
magnitude and the detection completeness for galaxies 
(see §4 and §5.1). After we placed the artificial galax- 
ies at random positions in the SSDF image with varying 
their fluxes, we measured their magnitude and isophotal 
area with the same technique as used for the observed 
galaxies. 

Figure lTEl shows a comparison between the model mag- 
nitude and the measured magnitude (total and isophotal) 
of the artificial galaxies. The measured total and isopho- 
tal magnitudes are almost consistent with the model 
magnitude within the uncertainties, although the mea- 
sured magnitude of the artificial galaxies at the faint end 
tend to be slightly fainter than the model magnitude. 
Thus, the extended wings of the observed PSF should 
not affect the measurement of the magnitude of the typ- 



FlG. 16. — Comparison between the model magnitude and the 
measured magnitude of the artificial gala xies that have the typical 
size in the SSDF image (see Figure [T^ . Filled circles show the 
total magnitude and filled triangles show the isophotal magnitude. 
Solid line shows where the model magnitude is equal to the mea- 
sured magnitude. The measured total and isophotal magnitudes 
are almost equal to the model magnitude. 

ical galaxies in the SSDF image. 

Figure llTf a) shows the isophotal area of the artificial 
galaxies with r^, = 0'.'25 as a function of the isopho- 
tal magnitude. We also show the isophotal area ver- 
sus isophotal magnitude diagram for small galaxies with 
Te — O'.'l in Figure llTr bl for a comparison. To esti- 
mate the effect of the extended wings of the observed 
PSF on the galaxy size measurement, we compared the 
isophotal area of the artificial galaxies convolved with the 
moffat profile and that convolved with the gaussian pro- 
file. Both of the profiles were fitted to the observed PSF. 
The moffat profile represents well the observed PSF that 
have extended wings, while the gaussian profile only rep- 
resents the central part of the observed PSF (see Figure 
114(1 . The isophotal areas of the = 0'.'25 galaxies con- 
volved with the moffat profile are almost consistent with 
those convolved with the gaussian profile at all magni- 
tude range. On the other hand, the isophotal areas of 
the Te = O'.'l galaxies convolved with the moffat profile 
are larger than that convolved with the gaussian pro- 
file. These implies that the extended wings of the ob- 
served PSF should not affect the isophotal area of the 
typical galaxies in the SSDF image, while it may affect 
the isophotal area of small galaxies at the bright magni- 
tude range. 

6.3. Effect of angular anisoplanatism 

The performance of AO correction degrades gradually 
with increasing distance from a guide star. This is due to 
increased decorrelation of the turbulence-induced aber- 
ration away from th e location of the guide star (angular 
anisoplanatism, e.g. lHardvll998D . Therefore, spatial res- 
olution (thus sensitivity) is usually highest in the region 
close to the guide star and decreases gradually as we go 
farther away from it. 

The angular anisoplanatism in our imaging data was 
evaluated from some snapshots of globular cluster M13 
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Fig. 17. — The isophotal area of the artificial galaxies are shown as a function of their isophotal magnitude (migo)- Left panel (a) shows 
the case for disk galaxies with effective radius (re) of 0'.'25, which is the typical size in the SSDF. We also show the same diagram for smaller 
disk galaxies with re = O'.'l in right panel (b) for a comparison. The artificial galaxies are convolved with a moffat profile (open circles) 
and a gaussian profile (crosses). The moffat profile represents well the observed PSF that have extended wings, while the gaussian profile 
only represents the central part of the observed PSF. The isophotal area of the re = 0'.'25 galaxies convolved with the moffat profile is not 
much different from that convolved with the gaussian profile. On the other hand, the isophotal area of the re = O'.'l galaxies convolved 
with moffat profile is larger than that convolved with the gaussian profile at rriiso < 23. 



using AO, in which a number of stars are distributed uni- 
formly in our 1' x 1' field of view. The snapshots were 
taken just after the deep imaging so that the observa- 
tional condition does not change significantly (Table 
The AO guide star for the snapshots was selected to have 
the same brightness as that of the deep imaging. Figure 
118b shows the variation of the stellar FWHM (spatial 
resolution) with respect to the distance from the guide 
star. Similarly, Figure IT^ shows the variation of the 5a 
limiting magnitude with an hour integration (sensitivity) 
in the field. The stellar FWHM and the limiting magni- 
tude deteriorate with increasing distance from the guide 
star. 

The stellar FWHM and the Scr limiting magnitude in 
the SSDF are O'.'IS and K' ~ 22.9 with an hour inte- 
gration, respectively, for the point-like source located at 
the distance of 24" from the guide star (dotted line in 
Figure [TH|l . Since more than 90% of the detected galax- 
ies in the SSDF are distributed in a range of distance 
from 10" to 35", the uncertainties in measurements of 
the stellar FWHM and the 5cr limiting magnitude due 
to anisoplanatism are about ±0'.'03 and ±0.1 mag, re- 
spectively, which should be negligible for the following 
discussions. 

7. DISCUSSION 

7.1. Galaxy number counts 

Figure 1191 shows our differential number counts of 
galaxies in the reliable magnitude range of 22 < K' < 25, 
where the signal-to-noise ratio for raw counts is greater 
than 3 and the detection completeness for point source 
is more than 50%. The count slope a = dlogN/dm is 
about 0.15 ± 0.08, where the uncertainty was defined as 
the uncertainty in the least-square fitting to the data. 
For a comparison, other K and K' band number counts 
in the literatures are also shown in Figure 1191 and sum- 
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Fig. 18. — Variation of (a) the FWHM of point sources (spatial 
resolution) and (b) the 5cr limiting magnitude with an hour inte- 
gration. The results as a function of distance from a guide star 
were based on the snapshots of a globular cluster M13 which were 
taken just after the deep imaging on 2003 April 23 (filled circles) 
and 24 (filled triangles) in Universal Time. Vertical dotted line 
shows the location of the reference point-like source ("S" in Figure 
from which the FWHM and 5cr limiting magnitude of the SSDF 
were derived (see §3). Most of detected galaxies in the SSDF are 
distributed in a range of distance from 10" to 35", and the effect 
of anisoplanatism should be negligible in this range. 



marized in Table|31 In this paper, we made no distinction 
between K and K' because the magnitude difference is 
expected to be {K' - K) < 0.03 (Minezaki et al. 1998a), 
which is negligible. 

Our counts at 23 < K' < 24.5 are in good a gree- 
ment with those of the SDF ijMaihara et al.ll200H) . al- 
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TABLE 3 

The comparison of the SSDF with other deep surveys in the K band. 



Survey 

This work (SSDF) 
Labbe et al. 2003 (FIRES) 
,Maihara ct al. 2001 (SDF) 
Bcrshady. ct al. 1998 
^ousteka^^L^^ 
^jorgovgk^^^^^995 



Instrument 


Integration time 


Area 






mag range"^ 




[hours] 


[arcmin^] 










Subaru/IRCS+AO 


26.8 


1.00 


25.0 


0.15 


22- 


-25 


VLT/ISAAC 


35.6 


6.25 


24.4 


0.15'i 


22- 


-24 


Subaru/CISCO 


9.7 


4.00 


24.4 


0.23 


20- 


-24 


Keck/NIRC 


4.9 


1.50 


24.0 


0.36 


19.5 


-22.5 


Keck/NIRC 


1.7 


2.41 


22.7 


0.23 


18- 


-23 


Keck/NIRC 


5.6 


0.44 


23.5 


0.32 


20- 


-24 



^Magnitude with 50% detection completeness for point source. 
^The slope of number counts a = dlog N / dm. 
*^Magnitude range for a. 

^Their slope is 0.25 at brighter magnitude of 20. 0< Ks <22.0 and declines to 0.15 at fainter magnitude 
range of 22. 0< Ks <24.0. 
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Fig. 19. — Differential X'-band number counts of galaxies es- 
timated from the SSDF galaxy counts in the reliable magnitude 
range (filled circles) where the signal-to-noise ratio for raw counts 
is more than 3 and the detection completeness for point source is 
more than 50%. Error bars for this work show the uncertainties 
due to a poisson statistics and a clustering of galaxies estimated 
from an angular correlation function. For a purpose of compari- 
son, the counts in the literatures are also plotted. Their magnitude 
ranges are limited by the same reliability criteria as ours. 



though these counts at 22 < K' < 23 differ by a factor 
of 1.2—1.6. Earher results of the if-band deep imag- 
ing by th e Keck telesc ope reported the slope of 0.23 
jMoustak as et all Il997fl 0.32 (Dioraovski ct al. 1995), 
and 0.36 ( Bershadv et al.ll99a|) at 20 < K' < 23, respec- 
tively. The discrepancy of the faint-end slopes could be 
originated from the difference in the adopted technique 
of incompleteness correction. The HDF-S deep imaging 
(FIRES, Labbe ct al. 2003) shows that the count slope 
changes from a ~ 0.25 at 20 < fsT' < 22 to a ~ 0.15 
at 22 < K' < 24. The break of the count slope at the 
faint end had not been seen in other deep i^-band imag- 
ing surveys. Although our counts are larger by a factor 
of about 1.5 than those of the HDF-S, our count slope 
at 22 < K' < 25 is consistent with the HDF-S slope at 
22 < K' < 24, supporting the flatter slope of a ~ 0.15 
at K' > 22. The break of the count slope at the faint 



end have also been seen in th e ff-band ga laxy counts in 
the northern HDF (HDF-N) l)ThomDSonl l20d3') that ex- 
tend to even fainter magnitude than ours. Assuming the 
average H — K color of the galaxies equals to 1, the H- 
band counts in the HDF-N show the break at the same 
magnitude as the if-band counts. 

Some theoretical models predicted further increase 
of number counts beyo nd K ~ 23. For example, 
iBabul fc FergusonI 1)19961) claimed that a large number 
of faint blue dwarf galaxies forming at high redshift 
shoul d beco m e obs ervable at this if-magnitude. More- 
over, iTomita ( 1995) claimed that the cosmic density at 
high redshift should be larger than that in the nearby 
universe, suggesting an excessive population of distant 
galaxies in the faint end. However, our number count 
slope is as flat as a ~ 0.15 down to K' ~ 25, suggest- 
ing that such extreme scenarios of high-redshift galaxy 
formation are unlikely. 

7.2. Extragalactic background light 

Extragalactic background light (EBL) in the optical 
and NIR wavelengths is an indicator of the total lu- 
minosity in the universe, which is believed to be dom- 
inated hy the integration of all stellar light ijBond et alJ 
119861 lYoshii fc Takahar'alll988() . Thereby, it provides a 
quantitative estimate of the baryonic mass that is a fun- 
damental quantity for galaxy formation and cosmology. 
If all stellar light is emitted from galactic systems, the 
EBL can be resolved into discrete galaxies by deep imag- 
ing sur vey. The faint-end s lope a of the HDF galaxy 
counts l)Williams et aljri996() is flatter than the critical 
slope of 0.4, with which the contributed flux from galax- 
ies to the EBL is constant against magnitude. There- 
fore, the light from the faint-end galaxies does not sig- 
nificantly increase the EBL, suggesting that the galax- 
ies that contribute to the E BL have already been re - 
solved into discrete galaxies (jMadau fc Pozzettil 12000(1 . 
Although there is a considerable scatter in the faint-end 
counts, the situation is the same for the NIR bands. 
However, the measurements of diffuse EBL in the op- 
tical and NIR bands suggest that the diffuse EBL flux is 
consistently higher than the integrated flux of galaxies. 
For instance, the K'-hand integ rated flux of galaxies was 
estimated at 7-8 nW m~^sr~^ l|Madau fc Pozz ctti 2000: 
iMaihara et al.ll200lT) . while the j^-b and diffuse EBL was 
estimated at 20-29 nWm-^sr-i l|Goriian et all 120001: 
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IWrighU I2OOII: iMatsumoto et al] l2001t iCambresv et al.l 
I2OOH) . This discrepancy could be caused by various 
selection effects or incompleteness of source detection. 
iTotani et al.l l|2001bf) theoretic a.lly estimated the nu mber 
of missed galaxies in the SDF l|Maihara e t al. 2001) due 
to selection effects or incompleteness of the detection of 
galaxies. They reported that the JsT-band EBL from the 
detected galaxies is estimated at 7.8—10.2 nWm~^sr~^ 
and concluded that more than 90% of the galaxies con- 
tributed to the EBL has already been resolved into dis- 
crete galaxies and the flux from missed galaxies cannot 
reconcile the discrepancy between the integrated flux of 
galaxies and the diffuse EBL. However, since there could 
be some inevitable uncertainties in theories, even deeper 
observation is essential to reveal the origin of this dis- 
crepancy. 

In this context, we derived the integrated flux of galax- 
ies contributed to the EBL from our if -band number 
counts, which is the deepest and less affected by the in- 
completeness of the source detection. Figure 1201 shows 
the count slopes in the separate magnitude ranges. At 
brighter magnitude, we used the slope of a = 0.59 ± 0.01 
at K' < 17.5 and a = 0.29 ± 0.01 at 17.5 < K' < 22, 
which were derived from the previous data with the 
signal-to-noise ratio of more than 3 and with the de- 
tection completeness for point source of more than 50%. 
At fainter magnitude of 22 < K' < 25, we used our re- 
sult of a = 0.15 ± 0.08, which is most reliable among 
the existing faint-end data. Figure |^ shows the con- 
tributed flux of galaxies in each ^fC'-magnitude bin to the 
EBL. The total contribution was calculated to be 9.26 
± 0.14 nWm~^sr~^ using the slopes in the above three 
magnitude ranges, where the uncertainty in the flux was 
determined from the uncertainty in the slope at the faint 
end. Even flatter count slope than the previous results 
at the faint end suggests that the contributed flux of the 
faint-end galaxies to the EBL is less significant. If the 
count slope does not change at K' > 25, the total con- 
tribution is expected to be about 9.41 nWm^^sr^^, sug- 
gesting that more than 98% of the galaxies contributed 
to the EBL has already been resolved as discrete galaxies 
in our data at K' < 25. 

Based on our SSDF counts down to K' 25, which 
is 0.5 mag deeper than the SDF, we found that our es- 
timated K'-hand EBL from galaxies is almost consistent 
with the theoretical estimation of Totani et al. (2001b). 
We also found that our estimated /-C'-band EBL from 
galaxies is close to the 7J-band EBL from galaxies es- 
timated from the HDF-N imaging (6 —7 nWm~^sr~^, 
iThompson et al.l2001l:lThompsonl2003(l . which IS compa- 
rable or even deeper than our j^T^- band imaging . These 
results support the conclusion of iTotani et al.l l)2001bfl 
that the EBL from galaxies has been resolved almost 
completely into discrete sources of normal galaxies. 
Comparison of our estimated EBL from galaxies with 
the measurements of diffuse EBL in the K band (20—29 
nWm~^sr~^) shows that a population of galaxies ac- 
counts for only less than 50% of the diffuse EBL in the 
K band. Unless the diffuse EBL measurements are con- 
taminated by crucial systematic uncertainties, a popula- 
tion other than known galaxy populations at even fainter 
magnitude is necessary to explain such discrepancy. It 
has been pointed out that a hypothetical population 
might consist of population III stars at very high redshift, 
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Fig. 20. — Three power-law curves were fitted to the differential 
number counts of galaxies at 12 < K' < 25 based on the combined 
data of ours and those in the literatures. The count slope a = 
d log N/dm is 0.59 ± 0.01 at K' < 17.5, 0.29 ± 0.01 at 17.5 < K' < 
22, and 0.15 ± 0.08 at 22 < K' < 25. 




e 10- 



12 14 16 18 SO 22 24 26 
K' total magnitude 



Fig. 21. — Contributed flux of galaxies in each i^'-magnitude 
bin to the extragalactic background light (EBL). The data in the 
literatures are also shown using the same symbols a s in Figure I^Hl 
Solid lines correspond to the fitted curves in Figured 

provided that a large amount of baryonic matter compa- 
rable to stars in galaxies should be used to forin such 
exotic stars at extremely high rate (Santos et al."2002 
ISalvaterra fc Fe rrara. 2003; Coorav & Yo shida.2004,) . 

7.3. Size distribution 

Figure|221shows the size distribution of detected galax- 
ies in the isophotal if'-magnitude {rriiso) versus isophotal 
area (Aiso) diagram for the SSDF and the SDF. The ob- 
served mean size is shown by filled circles. The sizes of 
individual SDF galaxies (crosses in FigureE^) are highly 
biased due to the large PSF size especially at faint mag- 
nitude of ifj'so > 24. However, because such bias is very 
small for the SSDF data thanks to the high spatial reso- 
lution with AO, the size of galaxies can be distinguished 
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Fig. 22. — Size distribution of the detected galaxies in the isophotal area versus isophotal magnitude diagram for (a) the SSDF and (b) 
the SDF. The SDF data were taken from|Totani et al. (2001a). Crosses show the individual galaxies. Filled circles show the mean size in 
each one magnitude bin. The horizontal error bars show the bin size around the mean and the vertical error bars show the Icr dispersion of 
the isophotal area of the detected galaxies. The gray shaded region is where the isophotal area is less than the expected area of point source 
for each isophotal magnitude in the SSDF or the SDF. The dashed line shows the maximum isophotal area for each isophotal magnitude 
with the surface brightness (SB) threshold (23.64 mag/arcsec^ for the SSDF, 24.1 mag/arcsec^ for the SDF). Thanks to the high spatial 
resolution due to AO, we can resolve the galaxies having the area size of less than 0.1 arcsec^ scale, which could not be done with the SDF 
image. The dot-dashed line in (a) shows the constant SB of 22.2 mag/arcsec'^. The constant SB line matches the observed size distribution 
within the uncertainty, although the observed size of the faint-end galaxies might be biased by the PSF. This suggests that the SB is not 
much different between high-redshift (faint) and low-redshift (bright) galaxies. 



from that of PSF even in the faint end. The slope of 
fTiiso — Aiso relation in the SSDF seems to be constant at 
20 < K'^,^ < 26, yielding d log A,,o/dTO,,o = -0.32±0.12. 
This slope is almost consistent with that of constant sur- 
face brightness {dlogAiso/drjiiso = —0.4) within the un- 
certainties, while the slope in the SDF seems to be flatter 
in the faint end due to the large bias of the PSF. We note 
that the constant surface brightness of 22.2 mag/arcsec^ 
(dot-dashed line in Figurel^^) matches the observed size 
distribution with the confidence level of 95%. This im- 
plies that the surface brightness is not much different 
between high-redshift (faint) and low-redshift (bright) 
galaxies. 

Figure 123 compares the observed size distribution of 
galaxies to theoretical predictions from a standard model 
of pure luminosity evolution (PLE) of galaxies that as- 
sumes no number evolution or no merging of galaxies. 
We used the s ame m odel as Totani & Yoshii (2000) and 
Totani et all l)2001a(l . which is an updated version of 
^shuTifToQ^ . In this model, galaxies are divided into 
six morphological types (dE, gE, Sab, Sbc, Scd, Sdm) 
and their type-dependent luminosity evolutions are de- 
scribed by the model of Er™oto & Yoshii (1987) and 
lArimotoet al.l l|1993 l in which the star formation history 
is determined to reproduce the present-day colors and 
chemical properties of galaxies. All galaxies are simply 
assumed to be formed at a single redshift of zp=3. The 
surface brightness profiles of elliptical and spiral galax- 
ies are modeled by the de Vaucouleurs and exponential 
profiles, respectively. These profiles are convolved with a 
moffat profile which nicely represents the observed PSF 
(see §6.2), and are used to calculate the isophotal area. 
The effective radius of galaxies is determined from the 
observed luminosity (L) versus effective radius (rg) re- 
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Fig. 23. — Same figure in Figure l^^ al. but with prediction of 
the PLE model with no size evolution described in the text. Thin 
lines show the predicted mean size for six morphological types of 
galaxies (dE, gE, Sab, Sbc, Scd, Sdm). The thick solid line shows 
the predicted mean size for all types of galaxies. 

lation of local galaxies f Ben der et alJlT99^ 'Impev et aj 
1996). Detailed descriptio ns o f the m odel are gi ven in 
lYoshiii H199$) . iTotani fc Yoshiil l|2000|) and iTotani eT ai] 
(j2001aj) . The predicted rriiso — ^iso relations for different 
types of galaxies are shown in Figure EHl We found that 
the predicted mean size for all types (thick solid line) is 
in good agreement with the observed mean size (filled 
circles). 

Figure |^ compares the observed size distribution to 
a PLE model that allows for intrinsic size evolution, in 
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Fig. 24. — Same figure as Figure |22ja) , but with prediction of 
the PLE model which allows for size evolution, assuming that the 
effective radius of galaxies, re, evolves as re oc (1 + 2)£. SoUd line, 
which is identical to the thick solid line in Figure 1231 shows the 
predicted mean size for all types of galaxies assuming no intrinsic 
size evolution = 0). The long-dashed, dot-dashed, dotted, and 
short-dashed lines show the models with size evolution of C=l) 2, 
-1, and -2, respectively. 

which galaxy size changes without changing the total lu- 
minosity (i.e. a size evolution is not driven by a number 
evolution). We here introduce a phenomenological pa- 
rameter and multiply an additional factor (1 + z)^ to 
the size of PLE model galaxies used in Figure |^ It is 
seen from Figure that the observed size distribution 
agrees with the model of C = 0, which favors no intrin- 
sic size evolution, in agreement with the earlier claim by 
lYos hii (1993). This suggests that the L — relation of 
high-redshift galaxies is almost same as the local rela- 
tion. Therefore, the result in Figure that the surface 
brightness is not much different between high- and low- 
redshift galaxies is further strengthened. 

JTotani et al. ( 2001a) performed the same analysis as 
ours using the SDF data taken by the Subaru/CISCO. 
They claimed that the size distribution of galaxies has 
the same trend as ours, although it was difhcult to dis- 
tinguish between the models of different types or different 
C based on their faint-end data at Kl^^ > 24, because the 
size of detected galaxies in the SDF was almost close to 
the size of PSF with FWHM ~ 0'.'45. For our SSDF data 
with AO, the PSF size is small enough to be separated 
from the galaxies at K'^„ > 24, thereby confirming the 
claim by iTotani et al.l lj2001all unambiguously. 

We note that a size evolution driven by a number evo- 
lution or merging is not necessarily ruled out if such 
process proceeds without changing the magnitude ver- 
sus size relation in Figure |^ In fact, this relation is 
maintained if the typical surface brightness of galaxies 
is preserved by mergers. What we can conclude from 
our result is that the surface brightness of the faintest 
SSDF galaxies is not much different from that expected 
for a simple PLE model without size evolution. There- 
fore, our result places a strong constraint on the con- 
temporary CDM-based structure formation theory (e.g., 
iWhite fc Ree^ Il978t iNagashima et al. 2QQ2.) , in which 
galaxies are formed through merging of many smaller 



building blocks at different formation epochs. Further 
information can be obtained by the relation between red- 
shift and effective radius of galaxies, which should be 
obtained from our SSDF data through the galaxy pro- 
file fitting (see Figure El and the photometric redshift 
technique (Minowa et al. 2005, in preparation). 

8. SUMMARY 

We have presented the first result of deep ^fC'-band 
imaging of the SSDF with the Subaru AO. The achieved 
detection limit is K' ^ 24.7 {5a, 0'.'2 aperture) for point 
source and K' 23.5 (5(T, 0'.'6 aperture) for galaxies with 
the total integration time of 26.8 hours. The sensitivity 
gain due to the AO correction is about 0.5 mag. The 
achieved spatial resolution is O'.'IS in stellar FWHM on 
average. We found from the simulation that the extended 
wings of the observed AO PSF does not affect the mea- 
surement of the flux and size of the typical galaxies in the 
SSDF. Combination of high sensitivity and high spatial 
resolution by AO enables detailed morphological study of 
high-redshift galaxies with the profile fitting technique. 
From a similar AO observation of globular cluster M13, 
we confirmed that the variation of the spatial resolution 
and the sensitivity to the point source (anisoplanatism) 
is not significant within our I'x 1' field of view. 

The slope of differential galaxy number counts {a = 
dlogN/dm) is about 0.15 at 22 < K' < 25. Our 
count slope i s consistent with that of the HDF-S (FIRES, 
iLabbe et alJ |^03), in which the slope changes from 
a ~ 0.25 at 20 < if' < 22 to a ~ 0.15 at 22 < K' < 24. 
Thus, our data support the fiatter slope of a ~ 0.15 at 
K' > 22 with the deepest image to date, while the other 
X-band deep surveys show the steeper slope at the faint- 
end. The discrepancies of faint-end counts between ours 
and other surveys might be caused by either field-to-field 
variation or different technique of incompleteness correc- 
tion. 

The contributed ffux of galaxies to the EBL is about 
9.43 nWm~^sr~^, and 98% of which appears to be re- 
solved into discrete galaxies in the SSDF at K' < 25. 
However, our EBL estimate can account for only less 
than 50% of the total ffux of diffuse EBL which was 
independently measured by satellites in the similar K 
ba nd. Our resul t further strengthens the previous claim 
by iTotani et alJ (|2001b') that there must be an exotic 
population that significantly contributes to the EBL, un- 
less the measurements of diffuse EBL are contaminated 
by crucial systematic uncertainties. 

We examined the size distribution of galaxies in the 
isophotal magnitude versus isophotal area diagram. In 
this diagram, the size distribution of detected galax- 
ies was obtained down to the area size of less than 0.1 
arcsec^, which is less than a half of that of the previous 
SDF galaxies. We found that the relation of observed 
mean size as a function of magnitude is nicely explained 
by a simple PLE model. We examined a possibility of 
intrinsic size evolution using a PLE model that allows 
for size evolution, and found that the model with no size 
evolution gives the best fit to the data. We conclude 
that the surface brightness of galaxies at high redshift 
is not much different from that expected from the size- 
luminosity relation of present-day galaxies. We cannot 
rule out a size evolution associated with a number evo- 
lution of galaxies if their surface brightness is preserved. 
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In other words, our rcsuh of constant surface brightness 
should place a strong constraint on the merging process 
in a hierarchical structure formation theory. 
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